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The objectives of this investigation were to study the interaction of the atmosphere with 
the surface of Mars through the impact of small objects that would generate dust and set 
the dust into motion in the atmosphere. The approach involved numerical simulations of 
impacts and experiments under con trolled conditions. 

The principal results are: 

Numerical simulations of meteoroid entry into the atmosphere, their deceleration, 
ablation and disintegration due to aerodynamic, forces and subsequent impact onto the 
surface have been conducted for various meteoroid sizes, and velocities. 

1. Meteoroids with radii of about 5 m and larger decelerate little in the Martian 
atmosphere and hit the surface with typical energies for meteoroids with radii 5, 20 and 
100 m equivalent to the energy of 60 kt, 4 Mt and 500 Mt TNT respectively. The mass of 
the dust injected into the atmosphere exceeds the mass of the impactor by an order of 
magnitude. The dust cloud rises to high altitudes (as an example, it rises to 10 km, 20 km 
and 25 km at moments of 30 s, 20 s and 60 s respectively for a 100 m impactor). For all 
the impactors with radii less than 100 m the main energy remains in the atmosphere. 

2. Direct impact by large meteoroids on Mars is rare (impacts of 5-m, 20-m and 100 m 
meteoroids occurs on Mars every 2-4 years, once in about 300-800 years and once 5000 
- 20,000 years respectively). 

3. The velocities of the impulsive wind created by the impact-generated shock waves is 
well above the dust threshold. For example, the impulsive wind velocity is about 200 m/s 
at a distance of 2.5 km from the crater produced by a meteoroid with a radius of 20 m or 
at a distance of 0.5 km from the crater created by meteoroids with radius of 5 m. 

4. The fate of smaller meteoroids (radius of 3 m and less) substantially depends on their 
strength. Typically such meteoroids release energy in the air above the surface. For a 1 m 
meteoroid, the mass of dust lifted into the atmosphere is > 1 t, larger than the mass of a 
typical dust devils. Meteoroids with radii of 3 m impact Mars 1-2 times per year; radii of 
1 m - about 10-20 times a year. 

5. Meteoroids < 1 m mainly fragment in the atmosphere and are rather numerous (the 
number of impacts with sizes of 0.2 m is about 3000-10,000 per year). 

6. Numerical simulations enable the shape and size of the dust cloud to be modeled. The 
shape, size of plumes and fireballs change with the size of meteoroids, their velocity and 
angle of trajectory. For vertical impacts dust plumes resemble those of dust devils 
(mushroom-type clouds). 

7. A special 3D code was developed to describe the interaction of the wind with the dust 
plume. Results of simulations were compared with laboratory experiments. Numerical 
simulations demonstrate the formation of two antisymmetrical vortices formed in the air 
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around the plume. These vortices may be considered as the initial rotation for impact- 
generated dust devils. 

8. Detailed tables of thermodynamic and spectral properties of hot Martian atmosphere, 
meteoroid vapor, and Martian soil vapor and developed radiation hydrodynamic codes 
allow the determination of radiation signatures of impacts. For large impactors the 
radiation energy is mainly released in the phase of plume expansion from the ground. For 
small impactors the radiation energy is mainly released during entry through the 
atmosphere and its disruption due to aerodynamic forces. Thus, the shape of light curves 
of Martian bolides is similar to that of terrestrial bolides, but the altitude of the meteoroid 
fragmentation is lower than in the Earth atmosphere by about 30 km. 

9. Preliminary estimates show that electrification of the impact-generated dust cloud may 
influence the distribution of particles by sizes. Electric signatures created by impacts 
could help in detecting impacts from Martian rovers and orbiters. 

Details of some of these results are given in the attached publication: 

2002, Atmospheric disturbances and radiation impulses caused by large-meteoroid 
impacts in the surface of Mars, Salary System Res., 36(3), 175-192, Kosarev, I.B., 
T.V. Losseva, I.V. Nemtchinov, V.V. Shuvalov, and R. Greeley. 
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INTRODUCTION 

Impacts of planetesimals have played an important 
role in forming the atmospheres of terrestrial planets 
during the period of heavy bombardrnent of the mne^ 
Solar System (see, e.g., Ahrens et ai, 1989; Hunten 
etal 1989; Melosh, 1994), i.e., 0.5-1 billion years 
after the formation of planets of the Solar em^ 
Uter, impacts of asteroids and comets occurred with a 
lower frequency. However, even now they may have a 
Sreffect on the evolution of planetary atmo- 
spheres. 

At present. Mars is a cold and extreinely dry planet 
(Zurek 1992; Schofield et ai, 1997). However, space 
flights to Mars showed that at one time its atmosphere 
was orobablv much denser and wetter than it is in the 
current epoch (Pollack er a/.. 1987; 
et al 1992; Hunten, 1993; Nelson and Greeley, i999). 
Valleys similar to river valleys exist on Mars. Sand and 
dust might have been formed earlier as a result of ero- 
sion caused by moving water ^d glaciers or by the 
weathering of rocks similar to that observed on E^rth. 
Another piece of evidence for the past existence o 
water on Mars is the presence of clear-cut tehees of 
sedimentary deposits, which have been revealed in pho- 


tographs acquired by orbital stations with the help of 
high-resolution cameras (Malin and Edgett, 2WX)). 
These terraces are normally associated with the shores 
of lakes that have gradually dried up or shallovv seas. 
However, there are other explanations for the on gin ot 
these sedimentary deposits, namely, volcanic ash, or 
dust and fragments produced by meteoroid impacts, or 
a joint action of all these factors. The question of how 
much new dust can arise during meteoroid irnpacts 
remains to be clarified, and this is one of the problems 
we solve in this study. 

The study of the ejection and precipitation of dust 
and fragments after meteoroid impacts is of consider- 
able interest in this and other aspects, for example, in 
connection with the problem of the ongm of dust 
storms on Mars (Rybakov et ai, 1997). 

According to the laboratory experiments performed 
by Greeley et ai ( 1 980; 1 992), the typical wind velocity 
required to raise dust from a flat surface covered by 
20-600 pm dust grains is equal to 80-170 m/s. For the 
surface covered by agglutinate particles and small 
“pebbles," the characteristic wind velocity is 25-60 m/s 
at a pressure of about 5 mbars. 
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The characteristic wind velocities on the Martian sur- 
face at the Viking Lander-2 landing site during the period 
preceding the onset of the 1977 global dust storm did not 
exceed 8-16 m/s (Zurek, 1 992). The dust elevation mech- 
anism on Mars with such low wind velocities is the most 
incomprehensible aeolian process. 

The traditional explanation is to assume that the rise 
of dust is initiated by the eddy structures (“dust devils”) 
that originate in the unstable atmosphere (Greeley 
etai, 1981; Ryan and Lucich, 1983; Kahn et ai, 1992; 
Rennd et ai, 1998; Fernandez, 1998; Metzger et ai, 
1999; Barnes, 1999; Metzger, 1999; Muiphy, 1999). 
This rise is enhanced by the increase in air speeds on 
the periphery of “vertical eddies” (dust devils) and by 
the sucking effect at the centers of these eddies, which 
gradually move to a new place, thus involving new dust 
in the motion. Such structures have been revealed in the 
photographs taken at the Mars Pathfinder landing site. 
In one of the images, a column 80 m wide and 350 m in 
height, moving up along the slope of the Big Crater at 
a speed of 4.6-3.5 m/s, is clearly seen. The dust concen- 
tration (7.5 X lOy^ kg/m^) and the total mass (about 0.7 1) 
contained in this largest dust devil photographed on 
Mars have been estimated by measuring the opacity of 
this column. This is a rather small value compared to 
the mass of dust raised during the 1977 global dust 
storm (430 Mt) and the mass raised during a local storm 
( 1 3 Mt) (Martin, 1 995). We do not reject the possibility 
of the combined effect of a great number of dust devils, 
the theoretically justified possibility of forming very 
large eddies (Renno et ai, 1 998), and the role of some 
other factors (regional topographic features and so on). 
However, this mechanism is obviously unproved. It is 
partly for this reason that we addressed ourselves to the 
study of the dust elevation mechanism in impact pro- 
cesses as an alternative to the dust rise in convective 
vortex structures (dust devils). 

Some hold the view that the change in the density of 
the Martian atmosphere and its humidity occurred in a 
nonrnonotone way: there may have been returns to the 
conditions in which water appeared again and again on 
the surface or in the near-surface layer. Here, we will 
consider only impacts on the planet’s surface under 
atmospheric conditions close to the currently existing 
atmosphere, assuming that water is absent on the sur- 
face and in the near-surface layers. Analytical estimates 
and numerical calculations in which the sizes of mete- 
oroids are varied within certain limits also give an idea 
of the role of atmospheric density. 

There are many possible reasons for the change in 
the Martian climate (Kasting and Toon, 1989). Among 
them are the periodic variations in the obliquity of the 
Martian orbit with a characteristic scale on the order of 
10*-I0* years (Ward, 1973; Ward et ai, 1974) and the 
change in degassing intensity and volcanic activity 
(Mouginis-Mark etai, 1992; Owen, 1992). Finally, the 
Martian climate might have changed as a result of the 
decrease in the number of impacts of large meteoritic 


bodies with time. On the one hand, impacts of large 
meteoroids may cause the erosion of the atmosphere 
(Melosh and Vickery, 1 989; Vickery and Melosh, 1 990- 
Ahrens. 1993; Zahnle, 1993). On the other hand, they 
replenish the atmosphere by adding vapors from the 
meteoroid and from the material of the planet’s surface 
(Ahrens and O’Keefe, 1987; Melosh, 1989). 

According to Ahrens and O’Keefe (1987), Ahrens 
et ai, (1989), and Ahrens (1993), noticeable erosion of 
the planetary atmosphere begins at the “atmospheric 
breakthrough” due to the acceleration of the shock 
wave that penetrates into upper, more rarefied atmo- 
spheric layers (Kompaneets, I960; Andriankin et ai, 

1 962; Zel’dovich and Raizer, 1 966). For the Earth, such 
a breakthrough occurs for the impactor radius R = 
0.5 km and the impactor energy Ep = 4.5 x 10“ MtTfix 

The application of the model of the breakthrough of 
the exponential almo.sphere (Ahrens, 1993) to other 
planets gives the following scaling law: 

where H is the scale height of the exponential atmo- 
sphere, V„ is the escape velocity, and poo is the atmo- 
spheric density near the planet’s surface. The H values 
for the Earth and Mars are not very different (8 and 
1 1 km, respectively); the “runaway velocities" (escape 
velocities) differ to a greater extent (5 and 1 1 km/s). 
whereas the densities p^o differ by two orders of magni- 
tude. Hence, we find for Mars: /?p ~ 90 m and Ep ~ 250 Mt 
TOT. For bodies of this size and smaller, most of the 
impactor energy remains in the atmosphere, causing its 
strong disturbance. 

A considerable loss of the atmospheric mass due to 
its expulsion by the vapors of the meteoroid and surface 
layers of the planet occurs for Mars, according to 
Melosh and Vickery (1989) and Vicke^ and Melosh 
(1990), at much higher impactor energies, namely, at 
10* MtTNT, and for an impactor radius of -1 km. Here, 
we will consider bodies of considerably smaller size. 

Analytical or semianalytical models of the atmo- 
spheric breakthrough, which were mentioned above, as 
well as some numerical calculations (see, e.g., 
Andrushchenko et ai, 1981), are based on certain 
^sumptions that include, among others, the exponen- 
tial distribution of density with height and the point 
character of explosion, assumptions that are not justi- 
fied in reality. It is necessary, therefore, to resort to 
numerical calculations, using the real density distribu- 
tion and real laws of energy release in the atmosphere 
during the flight of the meteoroid and its impact on the 
surface. , 

The qualitative pattern of impact erosion of the 
atmosphere and the replenishment of the atmosphere 
by the vapors of the meteoroid and surface layers is 
confirmed by two-dimensional calculations made by 
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Shuvalov and Artemieva (2000). For the Earth, in the 
range of impactor sizes from 1 to 30 km and for impac- 
tor velocities from 20 to 50 km/s, the replenishment is 
quite substantial for both asteroids and comets, whereas 
atmospheric erosion is rather weak. However, the quan- 
titative results of numerical calculations differ appre- 
ciably from estimates based on simple models. For the 
Earth, at impact velocities of 20 km/s and for asteroid 
diameters of 30, 10, 3, and 1 km, the ratios of the mass 
loss from the atmosphere to the impactor mass amount 
to4.7xI0■^ 1.6 X 10^, 7.1x10^, and4.3xl0■^ respec- 
tively. A simple scaling recalculation to Mars leads to 
nearly the same relative losses by mass for impactors of 
3, 1, 0.3, and 0.1 km in radius. Certainly, the scaling 
provides only a rough estimate, showing, however, that 
atmospheric erosion for 1-m and lOO-m meteoroids 
considered here is insignificant. Most of the energy and 
mass remains in the atmosphere, resulting in strong, 
and sometimes large-scale and long-term, disturbances. 

At the present time, when the Martian atmosphere is 
rarefied, bodies of 1-5 m in size reach the surface of the 
“Red Planet,” in fact, without disruption (Rybakov 
etai, 1997; Nemtchinov et al., 1999b; Kosarev et ai, 
1999, 2(X)0). Such meteoroids with kinetic energy 
equal to an explosion energy on the order of 0.1-10 kt 
TNT disintegrate at altitudes of about 25-35 km 
(Tagliaferri et ai, 1994) with a rate of about 25 events 
per year (Nemtchinov et al . , 1 997). Taking into account 
the lower density of the Martian atmosphere, which 
corresponds, near the planet’s surface, to the density of 
the terrestrial atmosphere at an altitude of about 35 km, 
and the lower area of the Martian surface (which is 
smaller than the Earth’s surface area by about a factor 
of two), we obtain that the speed of impacts of such 
bodies on the Martian surface is five to ten events per 
year. Although the mean velocity of meteoroids impact- 
ing on the Martian surface is somewhat lower than the 
velocity of Earth-impacting meteoroids (Steel, 1985; 
Olsson-Steel, 1987; ^hnle, 1993), such meteoroids, or 
even slightly smaller ones, may cause strong local dust 
storms, which, under “favorable” conditions, may 
sometimes be transformed into regional or even global 
storms (Rybakov et al., 1997; Nemtchinov et al., 
1999b). 


THE MODEL USED FOR CALCULATING 
IMPACT PROCESSES. NUMERICAL MODELING 
OF GAS AND DUST MOTION 

The physical and gasdynamic processes that occur 
at impacts of meteoroids on the Martian surface are 
considered here by means of numerical modeling. Gas- 
dynamic calculations were carried out using the SOYA 
code (Shuvalov, 1999; Shuvalov et al., 1999), which 
enables us to simulate complex multidimensional 
streams with the special isolation of interfaces between 


different “materials” (soil, vapors, atmospheric gas, 
etc.) and to employ different types of equations of state 
(presented in analytical or tabular form). We also take 
into account the effect of gravity. 

In the calculations considered below, we used the 
ANEOS equations of state (Thompson and Lauson, 
1972; Melosh, 1989) for the condensed phase. The 
equations of state and optical properties of the Msutian 
atmosphere heated by a shock wave, as well as the 
equations and properties for the material of impact- 
evaporated near-surface layers of Martian soil and for 
the vapors ejected by meteoroids themselves, will be 
briefly outlined in the second part of this study. 

Conceptually, the SOVA method is similar to the 
CTH method used in the USA (McGlaun et ai, 1990). 
However we cannot make a detailed comparison in 
view of the absence of necessary information about the 
American program. 

In the model used in the present paper, we make pro- 
vision for the conversion of deformed solid soil into a 
set of isolated particles. The density of the material 
decreases during the motion of the substance around 
the crater (excavation stage). At the instant when the 
density in a certain computational cell drops below a 
certain critical value p^p the continuous substance is 
replaced by a set of discrete particles (dust grains, sand 
particles, “pebbles”) of identical mass. The initial 
velocity of fragments at the moment of their formation 
was taken to be equal to the velocity of the continuous 
medium at the instant of its fragmentation. Such a 
method has been used by Teterev (1999). Moreover, 
new aerosol particles may arise in the process of scat- 
tering and adiabatic cooling of the vapors of the impac- 
tor and surface layers due to condensation. These very 
small particles (1-10 pm in size according to 2^1dovich 
and Raizer, 1966) move at the same velocity as the sur- 
rounding gas and can therefore be described by the con- 
tinuous-medium model. 

The^ize distribution of particles ejected from the crater 
was fitted to the experimental data obtained at ground 
nuclear and TNT explosions on rocks (Adushkin and Spi- 
vak, 1993; Arkhipov et ai, 1997; Solov’ev and Shuvalov, 
1999). Apart from particles that were produced in impacts, 
fine dust grains (0.1-1 mm in size) that already existed on 
the Martian surface can be involved in the motion. An 
algorithm for calculating the elevation of dust from the 
surface will be described below. 

In order to analyze the subsequent evolution of the 
cloud of particles, we considered their motion with 
allowance for the heat and momentum exchange with 
the surrounding gas. We used the method of representa- 
tive markers, each describing the motion of a large 
number of particles (10^-10”, depending on their size) 
with equal velocities, temperatures, and trajectories 
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(Teterev, 1 999). The typical number of such representa- 
tives was 10^-10^. A special algorithm (Shuvalov, 
1999) provided the determination of trajectories and 
velocities for large fragments, which meet almost no 
resistance and move along purely ballistic trajectories, 
as well as for microparticles moving with the local 
velocity of the gas flow. This method also allows the 
determination of the dust sedimentation rate due to 
gravitation, which depends on particle size. Note that 
the characteristic time of relaxation (equalization) of 
the velocities of large fragments and microparticles, 
is given by 


where d is the size of particles, is the gas density, pj 
is the number density of particles, and u is the characteris- 
tic gas velocity. For the characteristic values u = 1 km/s, 
Ps= 2.5 g/cm’, and Pj = lO"^ glcm?, the times x, are 
equal to 0.0 1 , 0. 1 , 1 , and 1 0 s for = 0.0 1 , 0. 1 , 1 , and 
10 cm, respectively. 

It is evident that, for small dust grains (d = 0.01-1 cm), 
the time of velocity relaxation may be comparable to 
the time of the expansion of the plume for a meteoroid 
and crater of moderate size. In this case, a strong inter- 
action exists between dust and gas. Dust and particles 
move at nearly identical velocity. Near the crater, the 
mass of ejected material exceeds the air mass, and dust 
is hardly braked at all. As the plume size increases (with 
time or with an increase in the impactor size), the role 
of braking in the atmosphere is increased. However, for 
very large impactors, the atmospheric breakthrough 
occurs virtually without braking. 

In the case of a meteoroid impact on rock, the sizes 
of fragments may be large enough — several centime- 
ters or more for 1- to 10-m bodies. Such large frag- 
ments travel appreciable distances (several kilometers 
or even greater) without substantial braking. 

It should be kept in mind that, apart from the sub- 
stance ejected from the crater to fairly high altitudes in 
the form of separate fragments and dust, which then 
precipitate at sufficiently large distances from the 
impact point, so-called continuous ejecta mantles also 
arise. They are formed near the crater, predominantly in, 
the zone extending from the rim crest for distances on 
the order of one to two crater radii (Melosh, 1989). 
Since this substance moves at comparatively low speed; 
the atmospheric disturbances are also small, and we 
therefore do not consider these processes here and, 
even less so, the subsequent modification of the crater. 

The computational algorithm was tested by compar- 
ing results with experimental data, in particular with the 
data on a 500-t TNT ground explosion on Novaya Zem- 
lya (Solov’ev and Shuvalov, 1 999) and on some ground 
nuclear explosions (Adushkin and Spivak, 1993; 


Adushkin et ai, 1998), including the stage of dust pre- 
cipitation (Trubetskaya and Shuvalov, 1 994). 

RESULTS OF CALCULATIONS 
OF THE IMPACT OF A 5-M BODY 

For meter-sized and smaller bodies, ablation and 
fragmentation in the atmosphere, as well as the braking 
of these bodies and their fragments, are of great impor- 
tance. Larger bodies, e.g., 5-m-sized bodies, reach the 
Martian surface almost without braking and dismption 
(Rybakov et ai, 1997; Nemtchinov et ai, 1 999b; Arte- 
mieva and Shuvalov, 2001). The results of the numeri- 
cal modeling of the impacts of cosmic bodies were 
briefly described by Kosarev et al. ( 1 999) and Nemtch- 
inov etal. (20(X), 2001). These investigators considered 
the vertical impacts of bodies 5, 20, and 1 00 m in radius 
moving with the same velocity, namely, 20 km/s. Here, 
we discuss these results in greater detail. 

In the first case, it was assumed that the surface is 
covered by a thick layer of fine dust; however, the role 
of the porosity of the body, which facilitates evapora- 
tion, was ignored. The results of the calculation for this 
case showed that the soil ejection rate at the initial stage 
of crater formation (within 10-20 ms of the impact for 
a body 5 m in radius) is 1-3 km/s. The velocity of vapor 
corresponding to this stage is equal to 1 0 km/s. The cal- 
culations also showed that the “curtain” formed from 
ejected soil strongly affects the stream in accordance 
with estimates and laboratory experiments (Schultz and 
Gault, 1979; Schultz, 1992; Melosh, 1994; Bamouin- 
Jha and Schultz, 1998). 

The braking of a dense vapor in air is responsible for 
the development of Rayleigh-Taylor instability at the 
cloud boundary. This instability leads to the formation 
of large-scale eddies and to the mixing of dust and 
atmospheric gas. Within 200-300 ms of the impact, a 
jet of the impactor material carves through the target- 
vapor cloud. Later on, the entire mass of the meteoroid 
emerges from the crater. 

The results of the calculation of the distribution of 
the density of gas and dust grains at a time instant 3 s 
after the impact of a 5-m body are shown in Fig. 1 . Dots 
on the right are the markers that represent dust frag- 
ments. At this instant of time, the atmospheric inhomo- 
geneity scarcely affects the results of the calculation, 
because the shock wave and the cloud of dust and vapor 
rise only about 2.5 km above the surface. 

The parameters of the vapor cloud were calculated 
also for the impact of a body on a rocky surface. For 
large fragments, the ejection process resembles that of 
the case in which dust is completely ignored, because 
fragments interact only slightly with surrounding air. In 
the case of an impact on a thick dust layer, the plume 
geometry is largely determined by the moving dust. 
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Fig. 1. The distribution of the density of gas and dust particles 3 s after the impact of a meteoroid of radius 5 m at a speed of 20 km/s. 


The shape of the vapor cloud coincides with the shape 
of the dust “curtain.” 

Analysis of the dust cloud for both impact variants 
at later instants of time showed that the ultimate results 
are close to each other. At the later phase, the evolution 
of the cloud is largely determined by the Archimedean 
force and turbulent mixing. The dust cloud rises about 
10 km. Its horizontal size reaches several kilometers. 
Within 300 s of the impact, the mass of dust involved in 
motion is equal to 60 Mq» the impactor with mass Mq 
hits the dust layer, and equals 10 Mq in the case of an 
impact on a rocky surface. Thus, for a 5-m impactor 
(i.e., for an impactor with mass equal to 1.3 x 10^ t), the 
mass of injected dust comprises 1.3 x 10^-8 x ICf* t. In 
both cases, 5 min after the impact, the cloud is not 
transparent for incident solar radiation. 

Observations in the regions of the Viking and Mars 
Pathfinder landing sites did not show the presence of a 
thick dust layer. However, observations made from 
orbiters indicate that dust clustering (with a mean thick- 
ness of 1-6 m) takes place in some regions, for exam- 


ple, in dunes of up to 10-20 m in height, with a mean 
distance between neighboring dunes of 500 to 2400 m 
(Christensen and Moore, 1992; Greeley et al., 1992; 
Parker and Rice, 1997). In some craters, the dust layer 
also appears to be sufficiently thick. On average, how- 
ever, the dust layer is apparently not very thick. As to 
the stratified deposits of sedimentary rocks on Mars, 
which may have a significant thickness (that is, they 
may completely fill, up to the crests, rather large craters 
up to 200 km in size), their density, structure, and other 
properties are virtually unknown, except that these 
deposits are likely to be solidified (Malin and Edgett, 
2000; Edgett and Malin, 2(XX)). In view of this, we will 
henceforth use only the calculations of ejection corre- 
sponding to the impact on rocky material. The possible 
role of water and ice will be briefly discussed telow. 

Another factor that may influence the evolution of 
the dust cloud is the wake that follows the meteoroid 
during its flight through the atmosphere. Expanding 
vapor experiences lower resistance in a rarefied wake 
(Nemtchinov and Shuvalov, 1992); the latter can be 
considered as a kind of half-empty tube. In the case of 
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an impact on the dust layer, the stream proves to be 
essentially anisotropic if the wake is taken into account. 
In 300 ms, the shock wave induced by the explosion 
covers in the wake twice the distance that it does out- 
side the wake. At the 30-s point, the upper boundary of 
the shock wave reaches an altitude of - 1 5 km, and the 
dust cloud reaches an altitude of 10 km. Within 300 s of 
the impact, this cloud rises 15 km above the surface. 
However, only a small portion of the dust moves within 
the wake. As a result, the mass of dust in the plume 
turns out to be nearly the same as when there is no 
wake. 

IMPACT OF A 1-M METEOROID 

In the case of an impact on a rocky surface, the 
results of the calculation for a body of one size can be 
approximately recalculated for the body of another size 
on the basis of the geometric similarity: the sizes of the 
plume, crater (at the stage of its formation and excava- 
tion), fragments, shock wave, and wake, as well as all 
characteristic times, are proportional to the impactor 
size, whereas pressure, velocities, densities, and tem- 
peratures remain unchanged. This is approximately 
valid for moderate-sized bodies and up to certain 
instants of time* as long as the plume and the entire dis- 
turbed region develop in a uniform atmosphere (the 
characteristic dimensions of which are less than the 
scale height of the atmosphere), and provided that we 
neglect radiation effects (the role of radiation will be 
considered below). For example, the results of calcula- 
tions for /? = 5 m can be recalculated for the case of the 
impact of a body with /? = 1 m, if this body does not 
experience braking and disruption in the Martian atmo- 
sphere. The mass of injected, dust within 60 s of the 
impact should be as high as 1 00 1, which is considerably 
greater than the mass inferred from observations of typ- 
ical dust devils. 

The above considerations of geometric similarity 
were tested by direct calculations of the corresponding 
problem. The dust cloud reached an altitude of -2 km, 
and its horizontal size was estimated to be -1 km. By 
the 100-s point, the dust cloud contains 250 1 (or 25 Mq, 
where Mq is the impactor mass). Thus, similarity is only 
approximate. 

The calculations were also carried out for a lower 
impact velocity, namely, for 1 1 km/s. The cloud sizes 
were approximately the same, but the mass of dust was 
only ~80 1 , or 8 Mq. 

We have dealt so far with the impact of a body that 
was not disrupted in the atmosphere owing to its rather 
high strength. The ejection of particles from the grow- 
ing crater was assumed to be the major source of raised 
dust, which is also the case for contact nuclear and TNT 
explosions. However, in the case of contact explosions, 
a considerable portion of dust can be raised behind the 
shock wave (Glasstone and Dolan, 1977; Adushkin and 
Spivak, 1993). Simple estimates (Nemtchinov et al,. 


1999a, 1999b) show that in many cases relatively small 
bodies (of about 1 m in size) may actually be destroyed, 
braked, and may release energy mostly over the sur- 
face. We carried out calculations for the case where 
energy release occurs at a high rate and the develop- 
ment of the disturbed region resembles that of an explo- 
sion in air (Glasstone and Dolan, 1977) (Dust, which 
already existed on the surface, is blown away by air 
moving behind the shock wave, and the elevation of 
dust occurs in convective streams that originate during 
the rise of the hot region). These processes are similar 
in many respects to those that occur in the case of sur- 
face erosion by natural winds (Greeley and Iversen, 
1985; Greeley et al., 2000a). 

Comparison with nuclear explosions set off near the 
Earth’s surface (Glasstone and Dolan, 1977) shows that 
“wind erosion” dominates when the explosion is made 
at a relative altitude of 900 m/kt'^ above the Martian 
surface (we took into account the lower density of the 
atmosphere of Mars). 

Numerical Model of Rise of Dust from the Sur- 
face. The above-mentioned SOVA code (Shuvalov, 
1999) was used to model the rise of dust, which already 
exists on the surface, caused by impulsed winds. How- 
ever, in contrast to the previous calculations, we u.sed 
the Navier-Stokes equations instead of the Euler equa- 
tions. This is necessary in order to take into account tur- 
bulent viscosity and diffusion. We used the viscosity 
coefficient v and the diffusivity D in the following 
form: 

VXD = 0.2 

Here V and D are given in mVs and the flow velocity V, 
in m/s. This relationship allows us to reasonably 
describe the origination and evolution of the dust ped- 
estal in laige-scale explosions (Nemtchinov et al., 
1993; Adushkin and Nemtchinov, 1994). 

To calculate the diffusive spreading of dust, we 
employed a Monte Carlo method. The displacement br 
of a particle over the time interval t was determined as 
follows: 

5r=Vx + j{6Dx)''^, 

where j is the unit vector oriented in a random direction. 

The main problem with this model was the question 
of how to determine the flow of dust that rises from the 
surface to the atmosphfere. We consider an air layer with 
a thickness 5 = (DAr)''^, where At is the time step. On 
the one hand, the gas velocity in this layer is determined 
by the velocity of the free flow. On the other, the layer 
is braked as a result of the entrainment of particles by 
the gas flow. For time At, new dust particles diffuse into 
the layer of atmospheric gas with mass = 5po (per 
unit area) for distance 5. Braking disappears when the 
velocity drops below the critical velocity and when 
the blowing-away of particles stops. The momentum 
balance is determined by the equation 
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a5p,(V,- V„) = 

where pg is the air density, is the air velocity, a is the 
blowing-away efficiency, and AA/p is the mass of dust 
particles per unit area. Dust does not penetrate into the 
atmosphere if 

The diffusion and upward motion of gas particles 
are responsible for the spreading of the dust cloud to 
high altitudes. At the same time, a certain amount of 
dust settles down owing to the action of gravity, winds, 
and downward diffusion from dusty layers adjacent to 
the surface. However, the impact of particles on the sur- 
face leads to the rise of new dust. This “saltation” effect 
can be taken into account by adding the maiss PAA/^ to 
the mass AAfp where is the specific mass of particles 
falling on the surface AA/pj for time At, and P is the sal- 
tation efficiency. 

The proposed model contains several parameters: 
V„, a, and p. The velocity V„ depends on several fac- 
tors, including the surface structure, the size of dust 
grains, and roughness (Greeley et ai, 2000b). This 
quantity and two other constants (a, p) have been mea- 
sured in numerous laboratory experiments. However, 
these parameters are still unknown in the natural condi- 
tions on Mars. Here, in our preliminary studies, we var- 
ied these constants to evaluate their effect on the ulti- 
mate result. 

The Calculation of the Rise of Dust at the “Air 
Explosion’* Caused by the Disintegration and Rapid 
Braking of a 1-m Meteoroid. The model described 
above was applied to the evaluation of the amount of 
dust in the Martian atmosphere after disintegration and 
rapid braking of the meteoroid with an initial size of 
1 m and an initial velocity of 1 1 km/s at an altitude of 
100 m over the surface; braking leads to a short-term 
release of energy equivalent to the chemical energy of 
a 0.15 kt TNT explosion. 

The following parameters were used in one of the 
variants: a = 0.5, P = 0.5, and V„ = 50 m/s. The average 
size of dust grains was assumed to be 2 p.m, and their 
density was taken to be 2.5 g/cm^. A considerable por- 
tion of dust from the region within 1 km of the epicenter 
was involved in the motion. The height of the dust ped- 
estal was 5-10 m. A fireball, very similar to that pro- 
duced by a nuclear explosion, was formed (Glasstone 
and Dolan, 1977) and began to rise. It initially con- 
tained a negligible mass of dust. However, the vortex 
flow in the atmosphere, which originates during the rise 
of the fireball, causes the dust to move toward the epi- 
center and then upward. The dust column connecting 
the fireball with the dust pedestal is formed, and then 
dust from the column gradually fills up fireball, form- 
ing a toroidal cloud. During the first minute, the veloc- 
ity of the dust cloud was 20-30 m/s. The mass of dust 
raised over the first 2 s was 35 1 (or 3.5 Mg, where Mg is 
the meteoroid mass). By the 15-s point, the mass of dust 
contained in the toroidal cloud, the dust column, and 
the dust pedestal comprised, respectively, 1 6, 6, and 6 1 , 
and the total mass of dust reached 28 t (or 2.8 Mg). 


For comparison, we point out that an explosion of 
equivalent power on the Earth at the same relative alti- 
tude raises a cloud of dust with a mass of about 30 1. 

Erosion regions are smaller than the dust pedestal 
region, because the pedestal expands radially. Erosion 
reaches a maximum at a distance of KX)-2()0 m from 
the explosion point. 

We also undertook computations with other values 
of the parameters ot, P, and D. The mass of dust var- 
ies in inverse proportion to V„(in the range 20-50 m/s). 
The value reaches 95 t for V„ = 20 m/s. 

A decrease in the turbulent diffusion and viscosity 
coefficients by an order of magnitude lowers the mass 
of raised dust by only approximately two times. A two- 
fold change in the a coefficient and variations in the 
saltation coefficient P in the range 0-0.5 lead to a 
change in Af„ of 20-30% = 47 t for a = 1 and P = 

0.5, and Af„ = 23 1 when a = 0.5 and P = 0). An increase 
in the particle mass by an order of magnitude slightly 
(by 10-20%) changes the mass of raised dust. 

An explosive wave also arises at the impact of an 
undisrupted meteoroid on the Earth’s surface. In this 
case, detailed computations of erosion are still lacking, 
but the data presented above show that the mass of dust 
raised from the surface probably does not exceed the 
mass of dust ejected from the crater. 

Dust clouds of conical or funnel-like shape, which 
rise to an altitude of 1-6 km above the surface, were 
revealed in high-resolution photographs obtained by 
Viking Orbiter. They were interpreted as convective 
vortex structures — “dust devils” (Thomas and Gierash, 
1985; Greeley et ai, 1992). It should be noted that 
impact-induced dust structures have a similar shape and 
size. 

The process of transformation of convective plumes 
into vortex structures (dust devils) is yet to be under- 
stood. Kuzmin et ai (2(X)1) draws attention to the 
effects of the interaction of winds with obstacles of dif- 
ferent type (for example, with large boulders formed, in 
particular, at meteoroid impacts, with crater rims, etc.). 
Note that the dust column ejected on impact also pre- 
sents an obstacle for wind. The processes of interaction 
of this column with wind are substantially three-dimen- 
sional, and we contemplate investigating them in the 
future experimentally (in wind tunnels) and theoreti- 
cally (using the 3D version of the SOVA method). 

It is universally recognized that the dust storms on 
Mars develop under specific conditions in the atmo- 
sphere, which are close to the state of instability and 
favor dust motion as a result of the heating of the dust 
cloud by solar radiation. Recently, the possibility of 
developing a self-sustaining dust front caused by the 
collapse of a dust devil and the spreading of dense gas- 
dust jets over the surface was noted, especially in the 
presence of oblique relief (Parsons, 2000). A similar 
process should also take place at the collapse of an 
impact-induced plume. We plan to study these pro- 
cesses in a separate paper. 
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Fig. 2. The disirihution of the density of gas and dust panicles (left) and temperature (right) 0.? s after the impact of a meteoroid of 
radius 20 ni at a speed of 20 km/s. 


IMPACT OF A 20-M METEOROID 

We now consider the consequences of the impact of 
a large (20-m-sized) meteoroid on the surface at an 
impact velocity of 20 km/s. 

The distributions of gas density and dust grains at 
the 0.3-s point are presented in Fig. 2. The representa- 
tive markers on the left show the position of dust and 
fragments, as well as the distribution of air and vapor 
density. A large-scale turbulence caused by the devel- 
opment of Rayleigh-Taylor instability is clearly seen in 
this figure. 

Similar di.stributions within 3 s and 10 s of the 
impact are shown in Figs. 3 and 4 (left). The height of 
the rise of the shock wave at these times is still lower 
than the characteristic scale height of the atmosphere. 
We see that most of the dust (in the upper part of the 
plume) moves together with the shock wave. Strictly 
speaking, the shock wave is almost spherical in shape, 
but the sphere center is raised above the surface. This is 
related to the fact that the kinetic energy of the meteor- 
oid transforms initially into the energy of the com- 
pressed material of the cosmic body and the planet’s 
surface, then partly reverts to the kinetic energy of the 
jet, and only later transforms into the thermal and 


kinetic energy of the atmosphere and vapor cloud as a 
result of Jet braking. 

The shock wave travels more slowly along the sur- 
face than above it (the maximum cross section of the 
disturbed region in the horizontal plane is attained 
within 3 s of the impact at an altitude of about 2-3 km 
and within 10 s of the impact at an altitude of 4-6 km.) 
According to Figs. 3 and 4, dust and fragments begin to 
fall on the surface within the spot corresponding to the 
radius of the cloud of 2,5-3 km before the arrival of the 
shock wave. In our computations, we assume that dust 
falling on the surface is not reflected from it and does 
not entrain dust that is already present on the surface. In 
reality, saltation of new dust grains increases the dust 
concentration in the near-surface layer. The gas, which 
moves behind the shock wave along the surface with 
sufficiently high speed, may entrain dust from the sur- 
face layer. 

The position of the shock wave (rj and the gas 
velocity behind the wave on the surface are given in 
Table 1 . We see that the gas velocity is considerably 
higher than the critical velocity leading to the rise of 
dust in the rarefied Martian atmosphere at distances of 
tens of kilometers from the impact site. 
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The time dependence of the gas velocity at a dis- 
tance of 2.5 km from the impact site is given in Table 2. 
The velocity of impulsed wind exceeds the critical ele- 
vation velocity for at least 10 s. 

We chose the distance 2.5 km intentionally. This is 
the distance between the Mars Pathfinder landing site 
and the Big Crater with diameter of 1.0-1. 3 km, the 
greatest crater near the landing site, which is located to 
the southeast of this site. 

Aeolian effects, including wind erosion of rocks at 
Rock Garden, have been analyzed by Bridges et al. 
(1999), Greeley et al. (1999, 2000a), Basilevsky et al. 
(1999), Golombek et al. (1999), and Golombek and 
Bridges (2000). 

Rocks are largely eroded by southeasterly winds. At 
the same time, meteorologic data (Scchfield et al., 
1997; Murphey et al., 1998) and the data of computa- 
tions based on the global circulation model (Haberle 
et al., 1999) indicate that the speed of such winds at the 
present time is, as a rule, lower than 10 m/s, i.e., it is 
insufficient to lift particles. Strong winds blow predom- 
inantly from the northwest (Greeley et al., 1999). 
Hence, the conclusion was drawn that erosion is caused 
by paleowinds that once blew at a time of different cli- 
mate (Golombek and Bridges, 2000). We put forward 
the hypothesis that erosion of rocks might have been 
caused by impulsed winds induced by a meteoroid 
impact. Initially, it appeared natural that it was the 
impact that formed the Big Crater. 

However, traces of wind erosion caused by south- 
easterly winds were found not only on boulders at the 
Rock Garden, but also on the rims of 85% of the 
105 craters located within an area 4.75 x 6.04 km ne^ 
the Mars Pathjmder landing site, including the Big 
Crater itself (Kuzmin et al., 2001). Thus, it is appropri- 
ate to analyze the possibility of impulsed wind erosion 
due to the impact of a larger meteoroid that may have 
occurred somewhat farther away to the southeast. Anal- 
ysis of images of the Martian surface (similar to that 
performed by Kuzmin et al. (2001) but within a much 
larger area) that would indicate the possible crater-can- 
didate is yet to be done. Here, we turn to the determina- 
tion of parameters of the disturbed atmosphere and dust 
cloud for large meteoroids. 

IMPACTS OF LARGE METEOROIDS 

Let us consider the results of computations for a 
stony meteoroid with radius Rq — 100 m and an impact 
velocity of 20km/s. The impactor mass is equal to 1.3 x 
10’ t. The temperature and density distributions at the 
10-s point are presented in Fig. 5. Dots on the left-hand 
side of Fig. 5 show the distribution of particles of the 
ejected soil material. As can be se«n, the height of the 
conical plume at this time attains 20-25 km. The max- 
imum height of the shock wave and shock-compressed 
air (indicated in gray) at this time is 25 km. According 
to the equation of state and the Hugoniot relations, the 


“air” temperature in the shock wave is low at this stage, 
while vapors have already cooled down considerably as 
a result of adiabatic expansion, and their temperature is 
close to the phase transition temperature due to conden- 
sation, as with the case of an impact on the Moon — a 
body completely devoid of atmosphere (Nemtchinov 
etai, 1998a, 1998b). The effect of the braking and 
heating of vapors at the boundary of the shock-com- 
pressed layer hardly appears at all in this situation. 

Parameters of the shock wave, which travels along 
the surface (radiative heating is ignored) are listed in 
Table 3. On the left, we present the data obtained for a 
uniform atmosphere; on the right, data for a nonuni- 
form atmosphere. We used the values of density and 
temperature at different altitudes that were obtained 
during the descent of Mars Pathfinder (Magalhaes 
etai, 1999). As can be seen, in the case of a nonuni- 
form atmosphere, the shock-wave parameters on the 
surface are considerably lower than for a uniform atmo- 
sphere. 

The shape of the hot region at the 10-s point differs 
considerably from the shape obtained for the 20-m 
impactor at the 2-s point. These shapes would be iden- 
tical if the geometric-similarity law were valid; how- 
ever, in this case the similarity law does not hold, 
because the height of the plume exceeds the character- 
istic scale height of the atmosphere. 

The vapor cloud is elongated in shape. The slowing 
down of the plume containing the vapors of the shock- 
compressed layer and the related increase in tempera- 
ture in the shock wave are markedly weaker than for the 
impact of a 20-m body. The main part of the plume has 
a temperature of -0.2 eV, and a denser and hotter vapor 
jet exists only near the crater. The elongated shape of 
the plume, the increase in the cloud height, and the 
decrease in the cloud temperauire point to a consider- 
able violation of geometric similarity. (Otherwise, the 
fireball would be located somewhat closer to the sur- 
face, with the center at an altitude of about 10 km; how- 
ever, our calculations demonstrate that the plume center 
is located at an altitude of about 1 2 km.) 

The maximum radius of the shock wave at the 30-s 
point is -30 km, i.e., the average radial velocity is 
-1 km/s. As follows from Table 3, the radius of the 
shock wave on the Earth’s surface at this point is only 
15 km. However, the gas velocity behind the shock- 
wave front amounts to 50 m/s, i.e., exceeds the value 
sufficient to raise the dust. Thus, strong erosion might te 
expected within an extensive zone of at least 10-15 km in 
size during a time interval of no less than 20—30 s. We 
have not yet carried out detailed calculations of erosion 
for two reasons. First, the enormous sizes of the dis- 
turbed region and the constraints on the possible num- 
ber of computational points do not allow us to conduct 
the calculations with a sufficiently fine step for the 
entire computed region, including near-surface layers; 
these calculations require a modification of the compu- 
tational procedure, which is currently being worked 
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Fig. 3. The distribution of the density of gas and dust particles (left) and temperature (right) 3 s after tlie impact of a meteoroid of 
radius 20 m at a speed of 20 km/s. 
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Fig. 4. The distribution of the density of gas and dust particles iO s after the impact of a meteoroid of radius 20 m at a speed of 
20 km/s. 
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out. Second, radiation emitted from the hot region may 
cause additional effects that facilitate the rise of dust. 
These problems will be considered in the second part of 
this study. 

Finally, we should note that experiments are also 
necessary on the aeolian abrasion prt^uced by dust par- 
ticles with a higher speed and shorter duration of its 
action than those conducted by Kraft and Greely (2000). 

At the 60-s point (Fig. 6), the shock wave reached an 
altitude of -300 km, i.e., the mean velocity increased to 
5 km/s and exceeded the escape velocity for Mars. As a 
consequence, the effect of the atmospheric break- 
through begins to manifest itself, although only a small 
mass possesses such a high velocity. The cone of ejec- 
tion of the planet’s substance reaches an altitude of 
100-120 km, and the dispersion of particles occurs in 
an almost inertial way: the strong effect of gravity is yet 
not pronounced. 

By the 60-s point, the cone radius increased to -100 km 
(Fig. 6), which means that the maximum radial velocity 
of the cone even increased. As to the explosion prod- 
ucts. their maximum radius reaches 70 km, i.e., the 
radial velocity is less than 1.2 km/s at a maximum 
height of 1 20 km and a vertical velocity of -2 km/s. 

The limiting height of elevation of dust and soil 
fragments and their radial dispersion are estimated to 
be -1000 km, on the assumption that the ejecta move 
ballistically with the above-indicated velocities. The 
mass of the ejected material is equal to about 10 Mq, 
i.e., 10* t for the given meteoroid. Assuming that this 
material is uniformly distributed over the circle of the 
radius equal to this height, the average specific mass is 
equal to -30 g/m^ which corresponds to a layer thick- 
ness of -10 pm. 

The evolution of the dust cloud in the atmosphere at 
a later stage must be determined on the basis of global- 
circulation models. However, the impact of a large cos- 
mic body can appreciably modify the global circulation 
as a result of the absorption of solar radiation by ejected 
dust and the heating of the atmosphere caused by the 
incident plume and the damping of oscillations. 

IMPACT RATE 

An extensive program of telescopic observations 
provides information about the distribution of 1-km and 
larger cosmic bodies in circumterrestrial space. 
According to the estimates made by Steel (1985), the 
mean probability of collisions in a sample of 284 asteroids 
crossing the orbit of Mars is 3.5 x 10-'° per year. Since 
there are - 10^ such bodies with size > 1 km, the mean 
time interval between impacts is equal to 3 x 10^ years. For 
smaller bodies, extrapolation bas^ on the A/ ~ D-°- law, 
where a = 1.8— 2.0, and on lunar (Neukum and Ivanov, 
1994) and Martian (Strom et al., 1992; Ivanov, 1999) 
cratering data yields an impact probability equal to one 
impact event every 3000-5000 years for 100-m bodies 
and one event every 200-300 years for 20-m bodies. 


Table 1. The shock wave front position (R^) and the maxi- 
mum wind velocity (i/j) as a function of time for the impact 
of a 20-m meteoroid 


t, s 

Rs* kni 

Wj, m/s 

1 

0.70 

190 

2 

1.25 

150 

3 

1.70 

180 

4 

2.20 

210 

5 

2.65 

210 

6 

3.10 

210 

7 

3.60 

170 

10 

4.75 

150 

14 

6.15 

140 

20 

8.10 

no 


Table 2. Near-surface wind speeds behind the shock wave 
as a function of time al a distance of 2.5 km from the impact 
site of a meteoroid of radius 20 m 


r, s 

u, m/s 

4 

0 

4.5 

180 

5 

150 

6 

70 

7 

40 


Satellite observations of flashts on the entry of 
meteoroids into the Earth’s atmosphere provide satis- 
factory statistical data for bodies with sizes of up to 
1-3 m (Tagliaferri et ai, 1994; Nemtchinov et al., 
1997). These data can also be used for estimates of the 
rate of impacts on Mars, taking into account the differ- 
ences in the surface areas. Data of satellite observations 
for the Earth can be represented (Nemtchinov et al. 

1 997) in the form N* = 7£^ ’^ or iV* - where N* 
is the number of impacts per year and £ is the energy in 
kt TNT. A stony meteoroid of radius 5 m with an impact 
speed of 20 km/s has a kinetic energy of about 60 kt. 
This value slightly exceeds the maximum energy that 
was detected by satellites over several years of system- 
atic observations (40 kt TNT). The rate of impacts of 
such meteoroids on Earth, obtained by extrapolation of 
the above dependence, is one impact event every two to 
three years. The rate of impacts on Mars can be evalu- 
ated as half of the impact detection rate for the Earth 
(Adolfson etal., 1996). It amounts to one impact event 
every five to six years. 

Extrapolation to high 'energies is open to question, 
because it is in this size range that a kink in the IV*(D) 
dependence is possible. However, using this depen- 
dence for a 20-m body, we obtain an impact detection 
rate of one event every 300-400 years. Smaller bodies 
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Fig. 5. The distribution of the density of gas and dust particles (left)and temperature (right) 10 s after the impact of an asteroid 
of radius 100 m at a speed of 20 km/s. 


(e.g., those of 0.3-1 m in size) penetrate the Martian 
atmosphere more frequently, and for such bodies we 
should take into account fragmentation and braking in 
the atmosphere. However, even a low-altitude “explo- 
sion” may produce a local dust storm (Rybakov et ai, 
1997; Nemtchinov et aL, 1999b). For larger bodies 
(>1 km), such events occur more rarely and may have a 
strong effect on the atmospheric conditions during a 
long period of time. 

There are different estimates of the time of transi- 
tion to the “tenuous” and cold Martian atmosphere. The 


Table 3. Parameters of the shock wave propagating along the 
surface used in the nonuniform and uniform models of the 
atmosphere after an impact of a meteoroid of radius 100 m 


Uniform atmosphere 

Nonuniform atmosphere 

f, s 

R, km 

f/,km/s 

P>Po 

s 

/?, km 

U, km/s 

plPo 

2 

2.7 

0.85 

8.4 

3 

2.45 

0.13 

3.9 

6 

5.6 

0.56 

4.2 

10 

6.0 

0.09 

2.5 

20 

13.2 

0.33 

2.4 

20 

11 

0.06 

1.7 

60 

28 

0.16 

1.55 

30 

15 

0.05 

1.6 

. 200 

o\ 

oo 

0.11 

1.37 

54 

24 

0.04 

1.5 


upper estimate (3. 5-3. 8 Myr ago) corresponds to the 
end of the period of heavy bombardment. However, 
there are data indicating that a dense atmosphere and 
humid climate existed on Mars also after this period 
(Pollack etai, 1987), The lower estimate (50-100 Myr 
ago) is based on a hypothesis for recent volcanic activ- 
ity on Mars (Hartmann, 1998, 2000; Hartmann and 
Berman, 2000; Hartmann et ai, 1999). Nevertheless, 
meteoroid impacts on Mars, which possessed an atmo- 
sphere then similar to the current one, repeatedly 
occurred even during this relatively short period. 

During the period when Mars had a rarefied atmo- 
sphere (for instance, over a period of 1 0^-10® years, 
when the number of impacts of cosmic bodies with a 
radius of -100 m was 10^- 10"*), the overall mass of 
ejected fine dust grains and soil fragments reached 
10‘®-10‘^ t, which, with allowance for the Martian 
surface area (1.45 x 10^ km^) corresponds to the spe- 
cific mass of the global layer of sediments - 7 x 10“^- 
7 X 10~^ g/cm^. 

Precipitated du5t may be swept away into other 
regions of the planet, especially during the period of 
dust storms. According to laboratory experiments per- 
formed by Greeley et ai (2000b), up to 1.5 x 10“^g/cm^ 
may be carried away from a unit area per unit time. In 
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Fig. 6. The distribution of the density of gas and dust particles 60 s after the impact of an asteroid of radius 100 m at a speed of 
20 km/s. 


Other words, the layer with specific mass on the order 
of 10"^ g/cra^ can be transferred in one day, while the 
entire precipitated layer can be swept away in several 
days of a dust storm. Its place may be occupied by dust 
brought from other places. Dust precipitated on the sur- 
face is not only redistributed by winds over the surface 
but can be concentrated in certain places, for example, 
in craters and canyons. Moreover, chemical “cementa- 
tion” processes decrease the mobility of “old” dust. 
Regretfully, these processes are insufficiently studied to 
provide a reliable estimate of the characteristic time of 
dust immobilization. Nevertheless, our data demon- 
strate that meteoroid impacts can explain, generally 
speaking, the origin of the current global dust layer on 
Mars. 

Malin and Edgett (2000) are doubtful of the impact 
origin of thick layers of sediments filling some very 
extended regions of the Martian surface, because ejecta 
products must have been concentrated only near impact 
sites. Our calculations and estimates show that the zone 
of falling of dust and fragments after impacts may be 
sufficiently* wide. Moreover, an even greater surface 


mass can be reprocessed by rare impacts of larger mete- 
oroids (up to 1 km in size). We do not consider here 
impacts of such bodies, assuming that for them the 
spherical shape of the surface and atmospheric layers of 
the planet must be taken into account at a sufficiently 
early phase of evolution. 

ROLE OF WATER AND ICE 

When calculating impacts, we ignored the possibil- 
ity of the existence of inundated and ice structures at a 
depth below the surface of Mars, in the Martian mega- 
regolith (Squyres et al., 1992), whose mean porosity is 
equal to -25%. The characteristic depth of a porous 
layer on the Moon is 6.5 km. The corresponding depth 
for Mars is estimated to be 2.8 km. The depth of craters 
for the largest meteorqids considered here is on the 
order of, or less than, this value. 

Porosity near the surface of Mars may reach 
20-50% (Squires et al.,lS 92). Over at least the first bil- 
lion years, groundwater existed in the Martian crust, 
especially in its upper laver-.(megaregolith) disrupted 
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by impacts. At present, water apparently persists partly 
as ice in the cryolithosphere. On average, this amount 
of water is likely to match the water-layer thickness of 
about 100 m. Morphological features of the surface in 
some regions of Mars show that viscous creep of the 
Martian regolith is due to the deformation of ice in the 
soil. Based on the analysis of the structure of deposits 
around large craters, Squyres et al. (1992) arrived at the 
conclusion that ice-bearing layers lie at a depth of 
greater than 200-300 m and therefore do not affect the 
ejection of high-speed dust jets and fragments, which 
disturb the atmosphere at impacts of moderate-sized 
meteoroids. As to the impacts of sufficiently large 
meteoroids (100 m in size and larger), one should take 
into account the possible presence of groundwater and 
porosity of the surface layer in studying the processes 
of the formation and development of plumes and the 
aerodynamic processes caused by the origination of 
products of the disruption of the surface material in the 
atmosphere and by impact-induced atmospheric distur- 
bances. These factors not only facilitate the ejection of 
substance at impacts (see, e.g., Stewart et al., 2000), but 
also change the character of aeronomic processes as 
compared to the current processes (Barth et al., 1992) 
involving insignificant amounts of water. 

We come to the conclusion that analysis of distur- 
bances in the Martian atmosphere caused by meteoroid 
impacts is far from being completed. It is necessary to 
take into account in greater detail the surface erosion 
and the elevation of dust behind the shock wave, the 
formation of vortex structures due to the interaction of 
the impact-induced shock wave with the ballistic wave 
that originated during the motion through the atmo- 
sphere, the upward ejection of dust due to the “explo- 
sion of the surface layer” caused by the drop in pressure 
in the blast wave below the initial level, the interaction 
of impact-induced convective plumes with the dust col- 
umn arising after the impact, and some other factors. 
However, certain difficulties related to the lack of 
knowledge of the properties of the surface itself in var- 
ious regions arise in this context. The necessary data 
can be derived from observations of impacts occurring 
in the current epoch. In this case, we come up against 
an additional problem related to the properties of 
impacting bodies and their braking/disruption in the 
atmosphere. The detection of the radiation flash by 
orbital stations on the entry of the meteoroid into the 
atmosphere and after the impact, as well as the detec- 
tion of the rise of the dust cloud coupled with ground- 
based seismic data, might facilitate the search for fresh 
craters and aeolian effects around craters, as well as the 
analysis of images obtained. 

CONCLUSIONS 

1. Relatively small meteoroids reach the surface of 
Mars almost without disruption and braking and with 
much higher frequency than they reach the Earth’s sur- 
face. This is related to the very rarefied Martian atmo- 


sphere. The impact detection rate for -1-m bodies is 
several events per year. Bodies of 5, 20, and 100 m in 
size impact on the Martian surface with a frequency of 
one impact event every 5-6 years, every 300-800 years, 
and every 50(X)-20000 years, respectively. 

2. No atmospheric breakthrough occurs after the 
impacts of bodies of £ 100 m in size. Most of the 
energy and mass remains in the atmosphere, leading to 
strong, large-scale, and sometimes global and long- 
term, disturbances. 

3. The mass of dust entering the atmosphere at 
impacts of 1-m meteoroids exceeds the typical mass of 
a “dust devil.” Impacts of bodies of about 5, 20, and 
100 m in size give ri.se to the injection of dust with mass 
similar to that rai.sed during a local, regional, or even 
global dust storm. 

4. The main mechanism of dust entry into the atmo- 
sphere is the ejection of a newly formed dust from the 
growing crater and vapor condensation. However, the 
impulsed wind behind the front of the shock wave that 
travels along the surface leads to wind erosion and to 
the rise of the already existing dust with a mass compa- 
rable to the ma.ss of the new dust. 

5. It is shown that at considerable distances from the 
impact site (for a 100-m meteoroid, at di.stances up to 
10-15 km at least) the velocity of gas particles behind 
the shock wave front exceeds the velocity at which dust 
is effectively blown away. 

6. A hypothesis is put forward that the abrasive 
action on boulders located in the Rock Garden region 
and on rims around small craters near the Mars Path- 
finder landing site, as well as near the Big Crater with a 
diameter of -I km, is caused by impulsed wind behind 
the shock wave front, which blew from the site of 
impact of a large meteoroid located to the southeast of 
the landing site, rather than by north-westerly paleow- 
inds (in contrast to current northeasterly winds). 

7. Some factors may cause stronger erosion than 
that assumed in existing computational models. This 
necessitates a deeper insight into the relevant processes. 
The second part of this study is devoted to the determi- 
nation of the characteristics of the thermal radiation 
impulse and to the e.stimation of its effect on the rise of 
dust. 

8. The formation of dust at numerous impacts of 
large meteoroids, which occurred over a long period 
(10’- 10* years), may explain the thickness of the cur- 
rent global du.st layer. However, the global circulation 
of the atmosphere and the redistribution and immobili- 
zation of dust should be taken into account to provide a 
more exact answer to this question. 
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